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Cosmic strings, a hot subject in the 1980's and early 1990's, lost its appeal 
when it was found that it leads to inconsistencies in the power spectrum of the 
measured cosmic microwave background temperature anisotropics. However, 
topological defects in general, and cosmic strings in particular, are deeply 
rooted in the framework of grand unified theories. Indeed, it was shown that 
cosmic strings are expected to be generically formed within supersymmetric 
grand unified theories. This theoretical support gave a new boost to the field 
of cosmic strings, a boost which has been recently enhanced when it was 
shown that cosmic superstrings (fundamental or one-dimensional Dirichlet 
branes) can play the role of cosmic strings, in the framework of braneworld 
cosmologies. 

To build a cosmological scenario we employ high energy physics; inflation 
and cosmic strings then naturally appear. Confronting the predictions of 
the cosmological scenario against current astrophysical/cosmological data we 
impose constraints on its free parameters, obtaining information about the 
high energy physics we employed. 

This is a beautiful example of the rich and fruitful interplay between 
cosmology and high energy physics. 

1 Introduction 

The basic ingredient in cosmology is general relativity and the choice of a 
metric. The Friedmann-Lemaitre-Robertson- Walker (FLRW) model, known 
as the hot big bang model is a homogeneous and isotropic solution of Ein- 
stein's equations; the hyper-surfaces of constant time are homogeneous and 
isotropic, i.e., spaces of constant curvature. The hot big bang model is based 
on the FLRW metric 

ds^ = dt^ - a^{t)-i,jdx^dx^ = a'^(T)[dT^ - j.jdx'dx^] ; (1) 



a{t) or a^r) is the cosmic scale factor in terms of the cosmological time t or 
the conformal time r (with dt =^ adr) respectively, and ^ij is the metric of a 
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space with constant curvature k. The metric 7^ can be expressed as 

-fijdx'dx^ = dr^ + x^{r) {dd^ + sin^ddip^) , (2) 

{r for K = 

s'mr for k = 1 (3) 
sinh r for k — ~1 ; 

the scale factor a(r) has been rescaled so that the curvature is k = ±1 or 0. 

The cornerstone of the hot big bang model is the high degree of symme- 
try of the FLRW metric: there is only one dynamical variable, the cosmic 
scale factor. The FLRW model is so successful that became the standard 
cosmological model. The high degree of symmetry of the metric, originally a 
theorist's simplification, is now an evidence thanks to the remarkable unifor- 
mity of temperature of the Cosmic Microwave Background (CMB) measured 
first by the COBE-DMR satellite [T]. 

The four pillars upon which the success of standard hot big bang model 
lie are: (i) the expansion of the Universe, (ii) the origin of the cosmic back- 
ground radiation, (iii) the synthesis of light elements, and (iv) the formation 
of large-scale structures. However, there are questions, which mainly concern 
the initial conditions, to which the hot big bang model is unable to provide 
an answer. These shortcomings of the FLRW model are: (i) the horizon prob- 
lem, (ii) the flatness problem, (iii) the exotic relics, (iv) the origin of density 
fluctuations, (v) the cosmological constant, and (vi) the singularity problem. 
To address these issues, inflation was proposed |2|3j . Inflation essentially con- 
sists of a phase of accelerated expansion, corresponding to repulsive gravity 
and an equation of state 3p < —p, which took place at a very high energy 
scale. Even though inflation is at present the most appealing scenario to de- 
scribe the early stages of the Universe, the issue of how generic is the onset 
of inflation is still under discussion |4|5|6j . at least within a large class of 
inflationary potentials, in the context of classical general relativity and loop 
quantum cosmology. 

From the observational point of view, the remarkable uniformity of the 
CMB indicates that at the epoch of last scattering, approximately 2 x lO^yr 
after the big bang, when the Universe was at a temperature of approximately 
0.26eV ~ 3 X lO'^K, the Universe was to a high degree or precision (10^^) 
isotropic and homogeneous. At very large scales, much bigger than llOMpc « 
lO^^km, the Universe is smooth, while at small scales the Universe is very 
lumpy. The fractional overdensity at the time of decoupling between baryons 
and photons was 

(7)^ =^x(f)^^(10-^-10-'); (4) 

the constant C depends on the nature of density perturbations and it is 
C = 0(10 — 100). Then one asks the following question: how does a very 
smooth Universe at the time of decouphng became very lumpy today? 
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In the 1980's and 1990's, cosmologists had the following picture in mind: 
small, primeval density inhomogeneities grew via gravitational instability into 
the large inhomogeneities we observe today. To address the question of the 
origin of the initial density inhomogeneities, one needs to add more ingre- 
dients (namely scalar fields) to the cosmological model. This is where high 
energy physics enter the picture. Clearly, to build a detailed scenario of struc- 
ture formation one should know the initial conditions, i.e., the total amount 
of nonrelativistic matter, the composition of the Universe, the spectrum and 
type of primeval density perturbations. 

For almost two decades, two families of models have been considered chal- 
lengers for describing, within the framework of gravitational instability, the 
formation of large-scale structure in the Universe. Initial density perturba- 
tions can either be due to freezing in of quantum fluctuations of a scalar field 
during an inflationary period, or they may be seeded by a class of topological 
defects which could have formed naturally during a symmetry breaking 
phase transition in the early Universe. On the one hand, quantum fluctuations 
amplified during inflation produce adiabatic, or curvature fluctuations with 
a scale-invariant spectrum. It means that there are fluctuations in the local 
value of the spatial curvature, and that the fractional overdensity in Fourier 
space behaves as \6k\'^ oc k~^. If the quantum fluctuations of the inflaton field 
are in the vacuum state, then the statistics of the CMB is Gaussian j8ji9,. On 
the other hand, topological defects trigger isocurvature, or isothermal fluctu- 
ations, meaning that there are fluctuations in the form of the local equation of 
state, with nongaussian statistics and a scale-invariant spectrum. The CMB 
anisotropies provide a link between theoretical predictions and observational 
data, which may allow us to distinguish between inflationary models and 
topological defects scenarios, by purely linear analysis. The characteristics 
of the CMB anisotropy multipole moments (position, amplitude of acoustic 
peaks), and the statistical properties of the CMB are used to discriminate 
among models, and to constrain the parameters space. 

Many particle physics models of matter admit solutions which correspond 
to a class of topological defects, that are either stable or long-lived. Provided 
our understanding about unification of forces and the big bang cosmology 
are correct, it is natural to expect that such topological defects could have 
formed naturally during phase transitions followed by spontaneously broken 
symmetries, in the early stages of the evolution of the Universe. Certain types 
of topological defects (local monopoles and local domain walls) lead to disas- 
trous consequences for cosmology and thus, they are undesired, while others 
may play a useful role. We consider gauge theories, thus we are only inter- 
ested in comic strings, since on the one hand strings are not cosmologically 
dangerous (monopoles and domain walls are), and on the other hand they 
can be useful in cosmology (textures decay too fast). 

Cosmic strings are linear topological defects, analogous to flux tubes in 
type-II superconductors, or to vortex filaments in superfluid helium. In the 
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framework of Grand Unifies Theories (GUTs), cosmic strings might have been 

formed at a grand unification transition, or much later, at the clcctrowcak 
transition, or even at an intermediate one. These objects carry a lot of energy 
and they could play a role in cosmology and/or astrophysics. In the simplest 
case, the linear mass density of a cosmic string, denoted by ^, is equal to 
the string tension. Thus, the characteristic speed of waves on the string is 
the speed of light. For strings produced at a phase transition characterised 
by temperature Tc, one expects roughly ^ ^ T^. The strength of the gravi- 
tational interaction of cosmic strings is given in terms of the dimensionless 
quantity ^ (Tc/Mpi)^, where G and Mp\ denote the Newton's constant 
and the Planck mass, respectively. For grand unification strings, the energy 
per unit length is ~ lO^^kg/m, or equivalently, G/z ~ 0(10"^ — 10~^). 

Topological defects (global or local) in general, and cosmic strings in par- 
ticular, are ruled out as the unique source of the measured CMB temperature 
anisotropics. Clearly, one should then address the following question: which 
are the implications for the high energy physics models upon which the cos- 
mological scenario is based? This leads to the following list of questions: (i) 
how generic is cosmic string formation? (ii) which is the role of cosmic strings, 
if any? and (iii) which is a natural inflationary scenario (inflation is still a 
paradigm in search of a model)? These questions will be addressed in what 
follows. We will see that cosmic strings are generically formed at the end of 
an inflationary era, within the framework of Supersymmetric Grand Unified 
Theories (SUSY GUTs). This implies that cosmic strings have to be included 
as a sub-dominant partner of inflation. We will thus consider mixed models, 
where both the inflaton fleld and cosmic strings contribute to the measured 
CMB temperature anisotropics. Comparing theoretical predictions against 
CMB data we will find the maximum allowed contribution of cosmic strings 
to the CMB measurements. We will then ask whether the free parameters of 
supersymmetric inflationary models can be adjusted so that the contribution 
of strings to the CMB is within the allowed window. 

Finally, the recent proposal that cosmic superstrings can be considered 
as cosmic string candidates opens new perspectives on the theoretical point 
of view. More precisely, in the framework of large extra dimensions, long 
superstrings may be stable and appear at the same energy scale as GUT 
scale cosmic strings. 

In what follows we first discuss, in Section 2, topological defects in GUTs. 
We classify topological defects and we give the criterion for their formation. 
We then briefiy discuss two simple models leading to the formation of global 
strings (vortices) and local (gauge) strings, namely the Goldstone and the 
abelian-Higgs model, respectively. Next, we present the Kibble and Zurek 
mechanisms of topological defect formation. We concentrate on local gauge 
strings (cosmic strings) and give the equations of motion for strings in the 
limit of zero thickness, moving in a curved spacetime. We subsequently dis- 
cuss the evolution of a cosmic string network; the results are based on heavy 
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numerical simulations. We then briefly present string statistical mechanics 
and the Hagcdorn phase transition. We end this Section by addressing the 
question of whether cosmic strings are expected to be generically formed af- 
ter an inflationary era, in the context of supersymmetric GUTs. In Section 
3 we discuss the most powerful tool to test cosmological predictions of theo- 
retical models, namely the spectrum of CMB temperature anisotropics. We 
then analyse the predictions of models where the initial fluctuations leading 
to structure formation and the induced CMB anisotropics were triggered by 
topological defects. In Section 4, we study inflationary models in the frame- 
work of supersymmetry and supergravity. In Section 5 we address the issue of 
cosmic superstrings as cosmic strings candidates, in the context of braneworld 
cosmologies. We round up with our conclusions in Section 6. 

2 Topological Defects 

2.1 Topological Defects in GUTs 

The Universe has steadily cooled down since the Planck time, leading to a 
series of Spontaneously Symmetry Breaking (SSB), which may lead to the 
creation of topological defects [7], false vacuum remnants, such as domain 
walls, cosmic strings, monopoles, or textures, via the Kibble mechanism |10] . 

The formation or not of topological defects during phase transitions, fol- 
lowed by SSB, and the determination of the type of the defects, depend on 
the topology of the vacuum manifold A^„. The properties of A^^ are usually 
described by the fc**^ homotopy group nk{Mn), which classifies distinct map- 
pings from the fc-dimensional sphere S'^ into the manifold A^„. To illustrate 
that, let us consider the symmetry breaking of a group G down to a subgroup 
H of G . If Ain = G/H has disconnected components, or equivalently if the 
order k of the nontrivial homotopy group is = 0, then two-dimensional 
defects, called domain walls, form. The spacetime dimension d of the de- 
fects is given in terms of the order of the nontrivial homotopy group by 
d = 4 — 1 — k. If A4„ is not simply connected, in other words if A4n contains 
loops which cannot be continuously shrunk into a point, then cosmic strings 
form. A necessary, but not sufficient, condition for the existence of stable 
strings is that the first homotopy group (the fundamental group) 7ri(A4„) of 
Mm is nontrivial, or multiply connected. Cosmic strings are line-like defects, 
d = 2. If A4n contains unshrinkable surfaces, then monopoles form, for which 
k = 1, d = 1. li Mn contains noncontractible three-spheres, then event-like 
defects, textures, form for which fc = 3, d = 0. 

Depending on whether the symmetry is local (gauged) or global (rigid), 
topological defects are called local or global. The energy of local defects is 
strongly confined, while the gradient energy of global defects is spread out 
over the causal horizon at defect formation. Patterns of symmetry breaking 
which lead to the formation of local monopoles or local domain walls are ruled 
out, since they should soon dominate the energy density of the Universe and 
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close it, unless an inflationary era took place after their formation. Local 
textures are insignificant in cosmology since their relative contribution to the 
energy density of the Universe decreases rapidly with time . 

Even if the nontrivial topology required for the existence of a defect is 
absent in a field theory, it may still be possible to have defect-like solutions. 
Defects may be embedded in such topologically trivial field theories [T^. 
While stability of topological defects is guaranteed by topology, embedded 
defects are in general unstable under small perturbations. 



2.2 Spontaneous Symmetry Breaking 

The concept of Spontaneous Symmetry Breaking has its origin in condensed 
matter physics. In field theory, the role of the order parameter is played 
by scalar fields, the Higgs fields. The symmetry is said to be spontaneously 
broken if the ground state is characterised by a nonzero expectation value of 
the Higgs field and does not exhibit the full symmetry of the Hamiltonian. 



The Goldstone Model To illustrate the idea of SSB we consider the simple 
Goldstone model. Let be a complex scalar field with classical Lagrangian 
density 

c = - , (5) 

and potential V{(j)): 

F(0) = jA[00-r;2]2 , (6) 

with positive constants A, rj. This potential, Eq. ([6|), has the symmetry break- 
ing Mexican hat shape. The Goldstone model is invariant under the U(l) 
group of global phase transformations, 

^{x) e'^^ix) , (7) 

where a is a constant, i.e., independent of spacetime. The minima of the 
potential, Eq. ([6]), lie on a circle with fixed radius |0| — f]; the ground state 
of the theory is characterised by a nonzero expectation value, given by 

(0|<^|0) = rje^" , (8) 

where 6 is an arbitrary phase. The phase transformation, Eq. ([7|), leads to the 
change 9 —> + a, which implies that the vacuum state |0) is not invariant 
under the phase transformation, Eq. ([7]); the symmetry is spontaneously bro- 
ken. The state of unbroken symmetry with (O|0|O) = is a local maximum 
of the Mexican hat potential, Eq. All broken symmetry vacua, each with 
a different value of the phase 9 are equivalent. Therefore, if we select the 
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vacuum with 6 = 0, the complex scalar field 4> can be written in terms of two 
real scalar fields, 02, with zero vacuum expectation values, as 

= ?/+^(0i+f^2) . (9) 

As a consequence, the Lagrangian density, Eq. ([5]), can be written as 

^ = liMl f + liM2 f - \W4>\ + Ant . (10) 

The last term, £int, is an interaction term which includes cubic and higher- 
order terms in the real scalar fields (j)i,(f)2- Clearly, 0i corresponds to a massive 
particle, with mass ^/Xr] > 0, while 02 corresponds to a massless scalar par- 
ticle, the Goldstone boson. The appearance of Goldstone bosons is a generic 
feature of models with spontaneously broken global symmetries. 

Going around a closed path L in physical space, the phase 6 of the Higgs 
field (f) develops a nontrivial winding, i.e., AO = 27r. This closed path can be 
shrunk continuously to a point, only if the field cf) is lifted to the top of its 
potential where it takes the value = 0. Within a closed path for which the 
total change of the phase of the Higgs field is 2tt, a string is trapped. A 
string must be either a closed loop or an infinitely long (no ends) string, since 
otherwise one could deform the closed path L and avoid to cross a string. 

We should note that we considered above a purely classical potential, 
Eq. (I6|), to determine the expectation value of the Higgs field 0. In a more re- 
alistic case however, the Higgs field is a quantum field which interacts with 
itself, as well as with other quantum fields. As a result the classical potential 
y(0) should be modified by radiative corrections, leading to an effective po- 
tential ycff(0)- There are models for which the radiative corrections can be 
neglected, while there are others for which they play an important role. 

The Goldtone model is an example of a second-order phase transition 
leading to the formation of global strings, vortices. 



The Abelian-Higgs Model We are interested in local (gauge) strings (cos- 
mic strings) , so let us consider the simplest gauge theory with a spontaneously 
broken symmetry. This is the abelian-Higgs model with Lagrangian density 

C = D^cj^V^^cj, - \f,,F^'' ~ F(0) , (11) 

where is a complex scalar field with potential V^(0), given by Eq. and 
Ffiu = d^Ai, — d^Af^ is the field strength tensor. The covariant derivative 
is defined by 2?^ = — ieA^, with e the gauge coupling constant and A^^ 
the gauge field. The abelian-Higgs model is invariant under the group U(l) 
of local gauge transformations 

0(x) ^ e^"("V(a;) ; A^{x) ^ A^{x) + -d^a{x) , (12) 
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where a{x) is a real single- valued function. 

The minima of the Mexican hat potential, Eq. ([6]), lie on a circle of fixed 
radius \4>\ = rj, implying that the symmetry is spontaneously broken and the 
complex scalar field (j) acquires a nonzero vacuum expectation value. Following 
the same approach as in the Goldstone model, we chose to represent (j) as 

= ^/ + ;^, (13) 

leading to the Lagrangian density 

^ = lid^M^ - ^fi'^l - \f^.F>^^ + \hPA^A>^ + Ant , (14) 

where the particle spectrum contains a scalar particle (Higgs boson) with 
mass rus = "v/Ar? and a vector field (gauge boson) with mass = \f2er\. The 
breaking of a gauge symmetry does not imply a massless Goldstone boson. 
The abelian-Higgs model is the simplest model which admits string solutions, 
the Nielsen- Olesen vortex lines. The width of the string is determined by the 
Compton wavelength of the Higgs and gauge bosons, which is ^ and 
~ m~^, respectively. 

In the Lorentz gauge, dfj^A^ = 0, the Higgs field (j) has the same form as 
in the case of a global string at large distances from the string core, i.e., 

« T^e"'' , (15) 

where the integer n denotes the string winding number. The gauge field 
asymptotically approaches 

A,,^ -d^\n<j> . (16) 

le 

The asymptotic forms for the Higgs and gauge fields, Eqs. and (|16p 
respectively, imply that far from the string core, we have 

V^q^^O , F^, «0. (17) 

As a consequence, far from the string core, the energy density vanishes ex- 
ponentially, while the total energy per unit length is finite. The string linear 
mass density /i is 

^ . (18) 

In the case of a global U(l) string there is no gauge field to compensate the 
variation of the phase at large distances from the string core, resulting to a 
linear mass density which diverges at long distances from the string. For a 
global U(l) string with winding number n = 1 one obtains 



f 


'ld(j)' 


Is 


rdd 



2iTrdr « 2nrf-\n(^ ) , (19) 
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where S stands for the width of the string core and i? is a cut-off radius at 
some large distance from the string, e.g. the curvature radius of the string, 
or the distance to the nearest string segment in the case of a string network. 
The logarithmic term in the expression for the string energy mass density 
per unit length leads to long-range interactions between global U(l) string 
segments, with a force ^ rf / R. 

The field equations arising from the Lagrangian density, Eq. pT|) . read 

[d^ - ze^^)(9^ - ieA^^)cj) H- ^c^{c^^ - if) = 

d^F^"" - 2elm [^{d'' - ieA'')<p\ = . (20) 

The equations of motion can be easily solved for the case of straight, static 
strings. 

The internal structure of the string is meaningless when we deal with 
scales much larger than the string width. Thus, for a straight string lying 
along the z-axis, the effective energy-momentum tensor is 

ft =M'5(2;)%)diag(l, 0,0,1) . (21) 



2.3 Thermal Phase Transitions and Defect Formation 

In analogy to condensed matter systems, a symmetry which is spontaneoulsy 
broken at low temperatures can be restored at higher temperatures. In field 
theories, the expectation value of the Higgs field (j) can be considered as a Bose 
condensate of Higgs particles. If the temperature T is nonzero, one should 
consider a thermal distribution of particles/antiparticles, in addition to the 
condensate. The equilibrium value of the Higgs field (p is obtained by min- 
imising the free energy F — E — TS. Only at high enough temperatures the 
free energy is effectively temperature-dependent, while at low temperatures 
the free energy is minimised by the ordered state of the minimum energy. 

Let us consider for example the Goldstone model, for which the high- 
temperature effective potential is 

Kff(0,T) = m\TMf = where m^{T) = ^{T^ - Gr^^) . (22) 

The effective mass-squared term (T) for the Higgs field in the symmetric 
state {(j)) = 0, vanishes at the critical temperature = VErj. The effective 
potential is calculated using perturbation theory and the leading contribution 
comes from one-loop Feynman diagrams. For a scalar theory, the main effect is 
a temperature-dependent quadratic contribution to the potential. Above the 
critical temperature, m?(T) is positive, implying that the effective potential 
gets minimised at = 0, resulting to a symmetry restoration. Below the 
critical temperature, rn?{T) is negative, implying that the Higgs field has a 
nonvanishing expectation value. 
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Even if there is symmetry restoration and the mean value {(p) of the 
Higgs field vanishes, the actual value of the field 4> fluctuates around the 
mean value, meaning that (j) any given point is nonzero. The thermal 
fluctuations have, to a leading approximation, a Gaussian distribution, thus 
they can be characterised by a two-point correlation function, which typically 
decays exponentially, with a decay rate characterised by the correlation length 
^. The consequence of this is that fluctuations at two points separated by a 
distance greater than the correlation length ^ are independent. 

Kibble [10] was first to estimate the initial density of topological defects 
formed after a phase transition followed by SSB in the context of cosmology. 
His criterion was based on the causality argument and the Ginzburg tempera- 
ture, Tg, defined as the temperature below which thermal fluctuations do not 
contain enough energy for regions of the field on the scale of the correlation 
length to overcome the potential energy barrier and restore the symmetry. 



AF is the difference in free energy density between the false and true vacua. 

According to the Kibble mechanism, the initial defect network is obtained 
by the equilibrium correlation length of the Higgs field at the Ginzburg tem- 
perature. Consequently, laboratory tests confirmed defect formation at the 
end of a symmetry breaking phase transition, but they disagree with defect 
density estimated by Kibble. More precisely, Zurek [13|14j argued that the 
the relaxation time f (T), which is the time it takes correlations to establish 
on the length scale has an important role in determining the initial 

defect density. 

Let us describe the freeze-out proposal suggested by Zurek to estimate 
the initial defect density. Above the critical temperature Tc, the field starts 
off in thermal equilibrium with a heat bath. Near the phase transition, the 
equilibrium correlation length diverges 



where v denotes the critical component. At the same time, the dynamics 
of the system becomes slower, and this can be expressed in terms of the 
equilibrium relaxation timescale of the field, which also diverges, but with a 
different exponent 



The values of the critical components /i, v depend on the theory under consid- 
eration. Assuming, for simplicity, that the temperature is decreasing linearly. 



e{TG)AF{TG) ~ Tg ; 



(23) 




(24) 




(25) 




(26) 
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where the quench timescale fq characterises the coohng rate. 

As the temperature decreases towards its critical temperature Tc, the 
correlation length ^ grows as 

m - (^) \ (27) 

but at the same time the dynamics of the system becomes slower, 

fit) ^ (^^) ' . (28) 

As the system approaches from above the critical temperature, there comes a 
time |<| during the quench when the equilibrium relaxation timescale equals 
the time that is left before the transition at the critical temperature, namely 

f{i) = \i\ . (29) 

After this time, the system can no longer adjust fast enough to the change of 
the temperature of the thermal bath, and falls out of equilibrium. At time t, 
the dynamics of the correlation length freezes. The correlation length cannot 
grow significantly after this time, and one can safely state that it freezes to its 
value at time i. Thus, according to Zurek's proposal the initial defect density 
is determined by the freeze-out scale [13114) 

i^ai)-f^^^^'^'^ ■ (30) 

We note that the above discussion is in the framework of second-order phase 
transitions. 

The above prediction, Eq. (j30p , has been tested experimentally in a variety 
of systems, as for example, in superfluid ''He |15ll6j and '^He |17ll8j . and in 
liquid crystals [19120121] . Apart the experimental support, the Kibble-Zurek 
picture is supported by numerical simulations 22 j23ii24j and calculations us- 
ing the methods of nonequilibrium quantum field theory [25| . 

2.4 Cosmic String Dynamics 

The world history of a string can be expressed by a two-dimensional surface 
in the four-dimensional spacetime, which is called the string worldsheet: 

= x^'{C) , a==0, 1 ; (31) 

the worldsheet coordinates are arbitrary parameters chosen so that 

is timelike and spacelike (= cr). 

The string equations of motion, in the limit of a zero thickness string, are 
derived from the Goto-Nambu effective action which, up to an overall factor, 
corresponds to the surface area swept out by the string in spacetime: 

So[x^] = -fiJ V^d\ , (32) 
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where 7 is the determinant of the two-dimensional worldsheet metric ^ab^ 
7 = det(7afc) = ]^e'"'e^'^ -fabled , lab = gf_,i,x'^^x'^f, . (33) 

If the string curvature is smaU but not neghgible, one may consider an expan- 
sion in powers of the curvature, leading to the following form for the string 
action up to second order 



S=- j d'CV^ip - fiiK^'K'' + I32R) , (34) 

where /3i,/32 are dimensionless numbers. The Ricci curvature scalar i? is a 
function of the extrinsic curvature tensor (with A = 1,2), 

R = K^'^'^K^f, ~ K^K^ ■ (35) 

— 7°''A'^{,. Finite corrections and their effects to the effective action have 
been studied by a number of authors [26]-|29|. 

By varying the action, Eq. ([32|) . with respect to a;^(C''), and using the 
relation ^7 = ii'^''d-fab, where lab is given by Eq. ([33b ). one gets the string 
equations o motion: 

<a^'^+^.^.7"''<a<6=0, (36) 

where F^^^ is the four-dimensional Christoffel symbol, 

C - la^'igr^.a + CJra.y " 9ua.j) , (37) 

and the covariant Laplacian is 



1 



-1 



daiV^r'x'^t,) ■ (38) 



One can derive the same string equations of motion by using Polyakov's form 
of the action [30| 

S[x^',hab] = / V^h'^'labd^C , (39) 

where hab is the internal metric with determinant h. 

Including a force of friction Ff^"^ due to the scattering of thermal particles 
off the string, the equation of motion reads [3T| 

-'^ + Fl^.r^x^.x^] = F'^(ui, T, a) . (40) 

The force of friction depends on the temperature of the surrounding matter 
T, the velocity of the fluid transverse to the world sheet u'^ = u'^ — x^^x'^'^Ucr, 
and the type of interaction between the particles and the string, which we 
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represent by cr. Cosmic strings of mass per unit length ^ would have formed 
at cosmological time 

to - (Gfiyhpi , (41) 

where tpi is the Planck time. Immediately after the phase transition the string 
dynamics would be dominated by friction [31j . until a time of order 

^ {G^iyhp, . (42) 

For cosmic strings formed at the grand unification scale, their mass per unit 
length is of order ~ 10~^ and friction is important only for a very short 
period of time. However, if strings have formed at a later phase transition, 
for example closer to the electroweak scale, their dynamics would be domi- 
nated by friction through most of the thermal history of the Universe. The 
evolution of cosmic strings taking into account the frictional force due to the 
surrounding radiation has been studied in Ref. [32] . 

The string energy-momentum tensor can be obtained by varying the ac- 
tion, Eq. (|32p . with respect to the metric g^^,, 

T^'V^ = -2-^ = / d\^j'^''x'^,x';,S<^^Hx^ - x^C)) ■ (43) 

For a straight cosmic string in a flat spacetime lying along the z-axis and 
choosing = t, C'^ = the above expression reduces to the one for the 
effective energy- momentum tensor, Eq. (|21|) . 



Cosmic Strings in Curved Spacetime The equations of motion for 
strings are most conveniently written in comoving coordinates, where the 
FLRW metric takes the form 

ds^ = a^{T)[dT^ - dr^] . (44) 

The comoving spatial coordinates of the string, x(r, a) , are written as a 
function of conformal time r and the length parameter cr. We have thus 
chosen the gauge condition = r. For a cosmic string moving in a FLRW 
Universe, the equations of motion, Eq. (|36|) . can be simplifled by also choosing 
the gauge in which the unphysical parallel components of the velocity vanish, 

X • x' = , (45) 

where overdots denote derivatives with respect to conformal time r and 
primes denote spatial derivatives with respect to a. 

In these coordinates, the Goto-Nambu action yields the following equa- 
tions of motion for a string moving in a FLRW matric: 
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The string energy per unit cr, in comoving units, is e = ■\/x'2/(l — x^), im- 
plying that the string energy is /ia J eda. Equation (|46)) leads to 

e a 

One usually fixes entirely the gauge by choosing a so that e — 1 initially. 

Cosmic Strings in Flat Spacetime In flat spacetime spacetime, the string 
equations of motion take the form 

daiV^l'^'x^l) = . (48) 



We impose the conformal gauge 

x-x' ^0 , x^ +x'^ ^0 , (49) 

where overdots denote derivatives with respect to (^^ and primes denote 
derivatives with respect to C^. In this gauge the string equations of motion 
is just a two-dimensional wave equation, 

X - x" = . (50) 

To fix entirely the gauge, we also impose 

i ^ = C" , (51) 

which allows us to write the string trajectory as the three dimensional vector 
x(i7, t), where = cr, the spacelike parameter along the string. This implies 
that the constraint equations, Eq. ([49]), and the string equations of motion, 
Eq. ([501), become 

i • x' = 

x2+x'2 = 1 

X - x" = . (52) 

The above equations imply that the string moves perpendicularly to itself 
with velocity x, that a is proportional to the string energy, and that the 
string acceleration in the string rest frame is inversely proportional to the 
local string curvature radius. A curved string segment tends to straighten 
itself, resulting to string oscillations. 

The general solution to the string equation of motion in flat spacetime, 
Eq. ([51;), is 



a((T -t)+ b((T + t) 



(53) 



where a((T — t) and h{a + t) are two continuous arbitrary functions which 
satisfy 

a'2 = b'2 = 1 . (54) 
Thus, cr is the length parameter along the three-dimensional curves a(CT), b(CT). 
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Fig. 1. Illustration of string intersections; (a) string-string intersection in one point, 
leading to the formation of two new long strings via exchange of partners; (b) 
string-string intersections in two points, leading to the formation of two new long 
strings via exchange of partners, and one closed loop; and (c) self-string intersections 
leading to the formation of one long string and a closed string loop 136! • 

Cosmic String Intercommutations The Goto-Nambu action describes to 
a good approximation cosmic string segments which are separated. However, 
it leaves unanswered the issue of what happens when strings cross. Numerical 
simulations have shown that the ends of strings exchange partners, intercom- 
mute, with probability equal to 1. These results have been confirmed for 
global [33], local [3T, and superconducting IIS] strings. 

String-string and self-string intersections leading to the formation of new 
long strings and loops are drawn in Fig. [TJ Clearly string intercommuta- 
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tions produce discontinuities in x and x' on the new string segments at the 
intersection point. These discontinuities, kinks, are composed of right- and 
left-moving pieces travelhng along the string at the speed of light. 

2.5 Cosmic String Evolution 

Early analytic work [37j identified the key property of scaling, where at least 
the basic properties of the string network can be characterised by a sin- 
gle length scale, roughly the persistence length or the interstring distance 
if, which grows with the horizon. This result was supported by subsequent 
numerical work [38j . However, further investigation revealed dynamical pro- 
cesses, including loop production, at scales much smaller than ^ |39|40j . 

The cosmic string network can be divided into long (infinite) strings and 
small loops. The energy density of long strings in the scaling regime is given 
by (in the radiation era) 

PL = kfit^^ , (55) 

where k is a numerical coefficient (k = 20 ± 10). The small loops, their size 
distribution, and the mechanism of their formation remained for years the 
least understood parts of the string evolution. 

Assuming that the long strings are characterised by a single length scale 
^{t), one gets 



Thus, the typical distance between the nearest string segments and the typ- 
ical curvature radius of the strings are both of the order of ^. Early numer- 
ical simulations have shown that indeed the typical curvature radius of long 
strings and the characteristic distance between the strings are both compa- 
rable to the evolution time t. Clearly, these results agree with the picture 
of the scale-invariant evolution of the string network and with the one-scale 
hypothesis. 

However, the numerical simulations have also shown |40l41j that small- 
scale processes (such as the production of small closed loops) play an essen- 
tial role in the energy balance of long strings. The existence of an important 
small scale in the problem was also indicated [1^ by an analysis of the string 
shapes. In response to these findings, a three-scale model was developed 
which describes the network in therms of three scales, namely the usual en- 
ergy density scale ^, a correlation length ^ along the string, and a scale ( 
relating to local structure on the string. The small-scale structure (wiggli- 
ness), which offers an explanation for the formation of the small sized loops, 
is basically developed through intersections of long string segments. It seemed 
likely from the three-scale model that ^ and ^ would scale, with ( growing 




(56) 
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Fig. 2. Snapshot of a networlc of long strings and closed loops in the matter era [49| . 

slowly, if at all, until gravitational radiation effects became important when 
C/'C ~ 10"'' j43l44j . Thus, according to the three-scale model, the small length 
scale may reach scaling only if one considers the gravitational back reaction 
effect. Aspects of the three-scale model have been checked [H] evolving a 
cosmic string network is Minkowski spacetime. However, it was found that 
loops are produced with tiny sizes, which led the authors to suggest [45] 
that the dominant mode of energy loss of a cosmic string network is parti- 
cle production and not gravitational radiation as the loops collapse almost 
immediately. One can find in the literature studies which support [46j this 
finding, and others which they do not [47].[48]. 

Very recently, numerical simulations of cosmic string evolution in a FLRW 
Universe (see, Fig.[2|), found evidence [49] of a scaling regime for the cosmic 
string loops in the radiation and matter dominated eras down to the hun- 
dredth of the horizon time. It is important to note that the scaling was found 
without considering any gravitational back reaction effect; it was just the re- 
sult of string intercmmuting mechanism. As it was reported in Ref. [49], the 
scaling regime of string loops appears after a transient relaxation era, driven 
by a transient overproduction of string loops with length close to the initial 
correlation length of the string network. Calculating the amount of energy 
momentum tensor lost from the string network, it was found [49] that a few 
percents of the total string energy density disappear in the very brief pro- 
cess of formation of numerically unresolved string loops during the very first 
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timesteps of the string evolution. Subsequently, two other studies support 
these findings [50].[5T]. 



2.6 String Thermodynamics 

It is one of the basic facts about string theory that the degeneracy of string 
states increases exponentially with energy, 

d{E) - ef^"^ . (57) 

A consequence of this is that there is a maximum temperature Tmax = 1//3h, 
the Hagedorn temperature [5 2^53,54) . In the microcanonical ensemble the 
description of this situation is as follows: Consider a system of closed string 
loops in a three-dimensional box. Intersecting strings intercommute, but oth- 
erwise they do not interact and are described by the Goto-Nambu equations 
of motion. The statistical properties of a system of strings in equilibrium are 
characterised by only one parameter, the energy density of strings, p, 

E , , 

P^JJ^ (58) 

where L denotes the size of the cubical box. The behaviour of the system 
depends on whether it is at low or high energy densities, and it undergoes a 
phase transition at a critical energy density, the Hagedorn energy density pH- 
Quantisation implies a lower cutoff for the size of the string loops, determined 
by the string tension [i. The lower cutoff on the loop size is roughly pT^I"^ ^ 
implying that the mass of the smallest string loops is toq ~ /i^^^- 

For a system of strings at the low energy density regime (p <C /oh), all 
strings are chopped down to the loops of the smallest size, while larger loops 
are exponentially suppressed. Thus, for small enough energy densities, the 
string equilibrium configuration is dominated by the massless modes in the 
quantum description. The energy distribution of loops, given by the number 
dn of loops with energies between E and E + dE per unit volume, is [53154155] 

dn cx e-°'^E^^/'^dE [p < pn) , (59) 

where a = (5/2mo) ln(/3H/p). 

However, as we increase the energy density, more and more oscillatory 
modes of strings get excited. In particular, if we reach a critical energy density, 
Ph, then long oscillatory string states begin to appear in the equilibrium state. 
The density at which this happens corresponds to the Hagedorn temperature. 
The Hagedorn energy density pjj: achieved when the separation between the 
smallest string loops is of the order of their sizes, is Pc ~ TOq. At the Hagedorn 
energy density the system undergoes a phase transition characterised by the 
appearance of infinitely long strings. 
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At the high energy density regime (p ^ pc), the energy distribution of 
string loops is |53l54l55j 

dn ArnJ^E'^/^dE {p > pn) , (60) 

where A is a numerical coefficient independent of mo and of p. Equation (|60p 
implies that the mean-square radius R of the closed string loops is 

R ~ m-'/'i?i/2 , (61) 

meaning that large string loops are random walks of step ~ ^tiq ^- Equations 
(|60l) and ([61]) imply 

dn = A'R-^dR {p > ph) , (62) 

where A' is a numerical constant. From Eq. (f62|) one concludes that the 
distribution of closed string loops is scale invariant, since it does not depend 
on the cutoff parameter mo- 

The total energy density in finite string loops is independent of p. Increas- 
ing the energy density p of the system of strings, the extra energy E — En, 
where i?H — PaL^, goes into the formation of infinitely long strings, implying 

p - Pinf = const (p > Ph) , (63) 

where pinf denotes the energy density in infinitely long strings. 

Clearly, the above analysis describes the behaviour of a system of strings 
of low or high energy densities, while there is no analytic description of the 
phase transition and of the intermediate densities around the critical one, 
p ~ Ph- An experimental approach to the problem has been proposed in 
Ref. and later extended in Ref. [57 . 

The equilibrium properties of a system of cosmic strings have been studied 
numerically in Ref. |56j . The strings are moving in a three-dimensional flat 
space and the initial string states are chosen to be a loop gas consisting of 
the smallest two-point loops with randomly assigned positions and velocities. 
This choice is made just because it offers an easily adjustable string energy 
density. Clearly, the equilibrium state is independent of the initial state. The 
simulations revealed a distinct change of behaviour at a critical energy density 
Ph — 0.0172 ± 0.002. For p < pn, there are no infinitely long strings, thus 
their energy density, pinf , is just zero. For p > pn, the energy density in finite 
strings is constant, equal to pn, while the extra energy goes to the infinitely 
long strings with energy density pinf = p — pn- Thus, Eqs. ([50]) and ([55]) are 
valid for all p > pn, although they were derived only in the limit p 3> pn- At 
the critical energy density, p = pn, the system of strings is scale-invariant. At 
bigger energy densities, p > pn, the energy distribution of closed string loops 
at different values of p were found |56j to be identical within statistical errors, 
and well defined by a line dn/dE oc E~^/^. Thus, for p > pn, the distribution 
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of finite strings is still scale-invariant, but in addition the system includes the 
infinitely long strings, which do not exhibit a scale-invariant distribution. The 
number distribution for infinitely long strings goes as dn/dE oc 1/E, which 
means that the total number of infinitely long strings is roughly log{E — En). 
So, typically the number of long strings grows very slowly with energy; for 
p > pn there are just a few infinitely long strings, which take up most of the 
energy of the system. 

The above numerical experiment has been extended |57j for strings moving 
in a higher dimensional box. The Hagedorn energy density was found for 
strings moving in boxes of dimensionality = 3, 4, 5 [57] : 

( 0.172 ± 0.002 for = 3 
PH = < 0.062 ± 0.001 for = 4 (64) 
[ 0.031 ± 0.001 for = 5 

Moreover, the size distribution of closed finite string loops at the high energy 
density regime was found to be independent of the particular value of p for 
a given dimensionality of the box de- The size distribution of finite closed 
string loops was found [STl to be well defined by a line 

where the space dimensionality da was taken equal to 3, 4, or 5 . The statis- 
tical errors indicated a slope equal to — (H-(iB/2)±0.2. Above the Hagedorn 
energy density the system is again characterised by a scale-invariant distribu- 
tion of finite closed string loops and a number of infinitely long strings with 
a distribution which is not scale invariant. 



2.7 Genericity of Cosmic Strings Formation within SUSY GUTs 

The Standard Model (SM), even though it has been tested to a very high 
precision, is incapable of explaining neutrino masses |58I59I60] . An extension 
of the SM gauge group can be realised within Supersymmetry (SUSY) . SUSY 
offers a solution to the gauge hierarchy problem, while in the supersymmetric 
extension of the standard model the gauge coupling constants of the strong, 
weak and electromagnetic interactions meet at a single point, Mqut — (2 — 
3) X 10^^ GeV. In addition, SUSY GUTs provide the scalar field which could 
drive infiation, explain the matter-antimatter asymmetry of the Universe, 
and propose a candidate, the lightest superparticle, for cold dark matter. We 
will address the question of whether cosmic string formation is generic, in the 
context of SUSY GUTs. 

Within SUSY GUTs there is a large number of SSB patterns leading from 
a large gauge group G to the SM gauge group Gsm = SU(3)cX SU(2)lX 
U(1)y- The study of the homotopy group of the false vacuum for each SSB 
scheme determines whether there is defect formation and identifies the type of 
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the formed defect. Clearly, if there is formation of domain walls or monopoles, 
one will have to place an era of supersymmetric hybrid inflation to dilute 
them. To consider a SSB scheme as a successful one, it should be able to 
explain the matter/anti-matter asymmetry of the Universe and to account 
for the proton lifetime measurements |58j . 

In what follows, we consider a mechanism of baryogenesis via leptoge- 
nesis, which can be thermal or nonthermal one. In the case of nonthermal 
leptogenesis, U(1)b-l (B and L, are the baryon and lepton numbers, respec- 
tively) is a sub-group of the GUT gauge group, Gqut, and B-L is broken at 
the end or after inflation. If one considers a mechanism of thermal leptogen- 
esis, B-L is broken independently of inflation. If leptogenesis is thermal and 
B-L is broken before the inflationary era, then one should check whether the 
temperature at which B-L is broken — this temperature defines the mass of 
the right-handed neutrinos — is smaller than the reheating temperature. To 
have a successful inflationary cosmology, the reheating temperature should 
be lower than the limit imposed by the gravitino. To ensure the stability 
of proton, the discrete symmetry Z2, which is contained in U(1)b-l, must 
be kept unbroken down to low energies. Thus, the successful SSB schemes 
should end at GsmX Z2. Taking all these considerations into account we will 
examine within all acceptable SSB patterns, how often cosmic strings form 
at the end of the inflationary era. 

To proceed, one has to flrst choose the large gauge group Gqut- In 
Ref. [61] this study has been done in detail for a large number of simple 
Lie groups. Considering GUTs based on simple gauge groups, the type of 
supersymmetric hybrid inflation will be of the F-type. The minimum rank of 
Ggut has to be at least equal to 4, to contain the Gsm as a subgroup. Then 
one has to study the possible embeddings of Gsm in Ggut so that there is 
an agreement with the SM phenomenology and especially with the hyper- 
charges of the known particles. Moreover, the large gauge group Gqut must 
include a complex representation, needed to describe the SM fermions, and it 
must be anomaly free. In principle, SU(n) may not be anomaly free. We thus 
assume that all SU(n) groups we consider have indeed a fermionic represen- 
tation which certifles that the model is anomaly free. We set as the upper 
bound on the rank r of the group, r < 8. Clearly, the choice of the maximum 
rank is in principle arbitrary. This choice could, in a sense, be motivated 
by the Horava-Witten [62] model, based on Eg x Eg. Concluding, the large 
gauge group Gqut could be one of the following: SO(IO), Eg, S0(14), SU(8), 
SU(9); flipped SU(5) and [SU(3)]^ are included within this list as subgroups 
of SO(IO) and Eg, respectively. 

A detailed study of all SSB schemes which bring us from Ggut down to 
the SM gauge group Gsm, by one or more intermediate steps, shows that 
cosmic strings are generically formed at the end of hybrid inflation. If the 
large gauge group Ggut is SO(IO) then cosmic strings formation is unavoid- 
able [SB. 
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The genericity of cosmic string formation for Eg depends whether one 
considers thermal or nonthermal leptogenesis. More precisely, under the as- 
sumption of nonthermal leptogenesis, cosmic string formation is unavoidable. 
Considering thermal leptogenesis, cosmic strings formation at the end of hy- 
brid inflation arises in 98% of the acceptable SSB schemes [63]. Finally, if 
the requirement of having Z2 unbroken down to low energies is relaxed and 
thermal leptogenesis is considered as the mechanism for baryogenesis, cosmic 
string formation accompanies hybrid inflation in 80% of the SSB schemes. 

The SSB schemes of SU(6) and SU(7) down to the Gsm which could ac- 
commodate an inflationary era with no defect (of any kind) at later times are 
inconsistent with proton lifetime measurements. Minimal SU(6) and SU(7) 
do not predict neutrino masses |61| . implying that these models are incom- 
patible with high energy physics phenomenology. 

Higher rank groups, namely S0(14), SU(8) and SU(9), should in general 
lead to cosmic string formation at the end of hybrid inflation. In all these 
schemes, cosmic string formation is sometimes accompanied by the formation 
of embedded strings. The strings which form at the end of hybrid inflation 
have a mass which is proportional to the inflationary scale. 



3 Cosmic Microwave Background Temperature 
Anisotropies 



The CMB temperature anisotropies offer a powerful test for theoretical mod- 
els aiming at describing the early Universe. The characteristics of the CMB 
multipole moments can be used to discriminate among theoretical models 
and to constrain the parameters space. 

The spherical harmonic expansion of the CMB temperature anisotropies, 
as a function of angular position, is given by 



dT 



/St 
df2n — {n)Y^Un) ; (66) 



yVe stands for the ^-dependent window function of the particular experiment. 
The angular power spectrum of CMB temperature anisotropies is expressed 
in terms of the dimensionless coefficients Ci, which appear in the expansion 
of the angular correlation function in terms of the Legendre polynomials Pf. 



ST ST 
^(n)^(n ) 



^ ^(2£ + l)QP,(cos^?)W2 ^ (g7) 



where we have used the addition theorem of spherical harmonics, i.e.. 
^ 21+1 



y,™(n)y/„(n') = -^P.(n • n' 



(68) 
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It compares points in the sky separated by an angle "(9. Here, the brackets 
denote spatial average, or expectation values if perturbations are quantised. 
Equation (I67|) holds only if the initial state for cosmological perturbations of 
quantum-mechanical origin is the vacuum [8|9| . The value of Ci is determined 
by fluctuations on angular scales of the order of 7r/£. The angular power 
spectrum of anisotropics observed today is usually given by the power per 
logarithmic interval in ^, plotting i{i + l)Ci versus 

To find the power spectrum induced by topological defects, one has to 
solve, in Fourier space for each given wave vector k a system of linear per- 
turbation equations with random sources: 

= 5 , (69) 

where T) denotes a time dependent linear differential operator, X is a vector 
which contains the various matter perturbation variables, and S is the ran- 
dom source term, consisting of linear combinations of the energy momentum 
tensor of the defect. For given initial conditions, Eq. ([55)1 can be solved by 
means of a Green's function, CJ(t, t'), in the form 

/To 
dr^jrn(To,r,k)5™(r,k) . (70) 

To compute power spectra or, more generally, quadratic expectation values 
of the form (Xj(ro, k)X^(To, k')), one has to calculate 

(X,(ro,k)X*(ro,k')) = 

/To To 
drg,„(T,k) j dT'^;;(r',k') X (5„(T,k)5:(T',k')) . (71) 

Thus, to compute power spectra, one should know the unequal time two-point 
correlators (5m (r, k)5*(T', k')) in Fourier space [64]. This object is calculated 
by means of heavy numerical simulations. 

The CMB temperature anisotropics provide a powerful tool to discrimi- 
nate among inflation and topological defects. On large angular scales {I < 50), 
both families of models lead to approximately scale-invariant spectra, with 
however a different prediction regarding the statistics of the induced per- 
turbations. Provided the quantum fields are initially placed in the vacuum, 
infiation predicts generically Gaussian fiuctuations, whereas in the case of 
topological defect models, the induced perturbations are clearly nongaus- 
sian, at least at sufficiently high angular resolution. This is an interesting 
fingerprint, even though difficult to test through the data. In the context of 
inflation, nongaussianity can however also be present, as for example in the 
case of stochastic inflation 65], or in a class of inflationary models involving 
two scalar flelds leading to nongaussian isothermal fluctuations with a blue 
spectrum [66 . In addition, allowing nonvacuum initial states for the cosmo- 
logical perturbations of quantum-mechanical origin, one generically obtains 
a non-Gaussian spectrum [8.9j , in the context of single-fleld inflation. 
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On intermediate and small angular scales however, the predictions of in- 
flation are quite different than those of topological defect models, due to the 
different nature of the induced perturbations. On the one hand, the inflation- 
ary fluctuations are coherent, in the sense that the perturbations are initially 
at the same phase and subsequently evolve linearly and independently of each 
other. The subsequent progressive phase shift between different modes pro- 
duces the acoustic peak structure. On the other hand, in topological defect 
models, fluctuations are constantly induced by the sources (defects). Since 
topological defects evolve in a nonlinear manner, and since the random initial 
conditions of the source term in the perturbation equations of a given scale 
leaks into other scales, perfect coherence is destroyed. The predictions of the 
defects models regarding the characteristics of the CMB spectrum are: 

• Global 0(4) textures lead to a position of the first acoustic peak at 
£ ~ 350 with an amplitude ~ 1.5 times higher than the Sachs- Wolfe 
plateau [57] . 

• Global 0(N) textures in the large N limit lead to a quite flat spectrum, 
with a slow decay after t ~ 100 [68]. Similar are the predictions of other 
global 0{N) defects [SHO]. 

• Local cosmic strings simulations [71] found a broad peak aX t k, 150 — 400, 
being produced from both vector and scalar modes, which peaks at ^ ~ 
180 and £ w 400 respcetively. 

The position and amplitude of the acoustic peaks, as found by the CMB 
measurements [72173174175] , are in disagreement with the predictions of topo- 
logical defect models. As a consequence, CMB measurements rule out pure 
topological defect models as the origin of initial density perturbations leading 
to the observed structure formation. 



3.1 Mixed Models 

Since cosmic strings are expected to be generically formed in the context 
of SUSY GUTs, one should consider mixed perturbation models where the 
dominant role is played by the inflaton field but cosmic strings have also a 
contribution, small but not negligible. Restricting ourselves to the angular 
power spectrum, we can remain in the linear regime. In this case, 

Ce = aCl + (1 - a)C| , (72) 

where and C| denote the (COBE normalized) Legendre coefficients due 
to adiabatic inflaton fluctuations and those stemming from the cosmic string 
network, respectively. The coefficient a in Eq. ()72|) is a free parameter giving 
the relative amplitude for the two contributions. Comparing the Ci, calcu- 
lated using Eq. (|72p — where is taken from a generic inflationary model 
and C| from numerical simulations of cosmic string networks — with data 
obtained from the most recent CMB measurements, one gets that a cosmic 
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Fig. 3. l)Ce versus £ for three different models. The string contribution turns 

out to be ~ 18% of the total [76] , 



string contribution to the primordial fluctuations higher than 14% is excluded 
up to 95% confidence level |76I77I78| (see, Fig.^. 

In what follows, we follow a conservative approach and do not allow cosmic 
strings to contribute more than 10% to the CMB temperature anisotropics. 

3.2 Supersymmetric Hybrid Inflation 

Inflation offers simple answers to the shortcomings of the standard hot big 
bang model. In addition, simple inflationary models offer successful candi- 
dates for the initial density fluctuations leading to the observed structure 
formation. 

One crucial question though is to answer how generic is the onset of 
inflation |4I5I6| and to find consistent and natural models of inflation from 
the point of view of particle physics. Even though one can argue that the 
initial conditions which favor inflationary models are the likely outcome of 
the quantum era before inflation [4j, one should then show that inflation 
will last long enough to solve the shortcomings of the standard hot big bang 
model |5I6| . In addition, to find natural ways to guarantee the fiatncss of the 
inflaton potential remains a difficult task. Inflation is, unfortunately, still a 
paradigm in search of a model. It is thus crucial to identify successful but 
natural inflationary models motivated from high energy physics. 

In what follows we discuss two well-studied inflationary models in the 
framework of supersymmetry, namely F/D-term inflation. Our aim is to check 
the compatibility of these models — here cosmic string inflation is generic — 
with the CMB and gravitino constraints. 
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F-term Inflation F-term inflation can be naturally accommodated in the 

framework of GUTs, when a GUT gauge group Ggut is broken down to the 
GsM at an energy Mqut, according to the scheme 

y^GVT > ill >rl2 ><JSM , (,'oj 

is a pair of GUT Higgs superfields in nontrivial complex conjugate 
representations, which lower the rank of the group by one unit when acquiring 
nonzero vacuum expectation value. The inflationary phase takes place at the 
beginning of the symmetry breaking Hi H2 . 

F-term inflation is based on the globally supersymmetric renormalisable 
superpotential 

W^^fl = kS{<1>+<1>. - M^) , (74) 

where 5* is a GUT gauge singlet left handed superfield, and are de- 
flned above; k and M are two constants (M has dimensions of mass) which 
can be taken positive with field redefinition. The chiral superfields S, <P+ , <P- 
are taken to have canonical kinetic terms. This superpotential is the most 
general one consistent with an R-symmetry under which W ^ e^^W ^ 
e~^^<l>_ ,<P^ e^^fp^^ and S e^^S. An R-symmetry can ensure that the 
rest of the renormalisable terms are either absent or irrelevant. 
The scalar potential reads 

V{ct>^A-,S) = \F^^\' + \F^_\' + \Fs\' + \ Y^glDl . (75) 

a 

The F-term is such that = \dW/d^i\0=o, where we take the scalar com- 
ponent of the superfields once we differentiate with respect to #i = ^-,S. 
The D-terms are 

Da = ^i{Tayjcly^ , (76) 

with a the label of the gauge group generators T^, ga the gauge coupling, 
and the Fayet-Iliopoulos term. By definition, in the F-term inflation the 
real constant is zero; it can only be nonzero if Ta generates an extra U(l) 
group. In the context of F-term hybrid inflation, the F-terms give rise to the 
inflationary potential energy density, while the D-terms are flat along the 
inflationary trajectory, thus one may neglect them during inflation. 

The potential has one valley of local minima, V = k^M^^ for S > M with 
0+ = = 0, and one global supersymmetric minimum, V — 0, ai S = Q and 
(j)^ — (f)_ = M. Imposing initially S » M, the fields quickly settle down the 
valley of local minima. Since in the slow roll inflationary valley, the ground 
state of the scalar potential is nonzero, SUSY is broken. In the tree level, 
along the inflationary valley the potential is constant, therefore perfectly 
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flat. A slope along the potential can be generated by including the one-loop 
radiative corrections. Thus, the scalar potential gets a little tilt which helps 
the inflaton field S to slowly roll down the valley of minima. The one-loop 
radiative corrections to the scalar potential along the inflationary valley, lead 
to an effective potential |79|80|81|82] 
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where yl is a renormalisation scale and J\f stands for the dimensionality of the 
representation to which the complex scalar components , of the chiral 
superfields 'P+,<P- belong. For example, Af = 27,126,351, correspond to 
realistic SSB schemes in SO(IO), or Eg models. 

Considering only large angular scales, i.e., taking only the Sachs- Wolfe 
contribution, one can get the contributions to the CMB temperature aniso- 
tropics analytically. The quadrupole anisotropy has one contribution coming 
from the inflaton field, splitted into scalar and tensor modes, and one con- 
tribution coming from the cosmic string network, given by numerical simu- 
lations [83. . The infiaton field contribution is 
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where the quadrupole anisotropy due to the scalar and tensor Sachs- Wolfe 
effect is 
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respectively, with V = dV{if)/dip, Mpi the reduced Planck mass, Mpi — 
(87rG)-i/2 ^ 2.43 x lO^^GeV, and ipq the value of the inflaton fleld when 
the comoving scale corresponds to the quadrupole anisotropy became bigger 
than the Hubble radius. It can be calculated using Eqs. (|77)) . (|78|l . ((79|) . 

Fixing the number of e-foldings to 60, the inflaton and cosmic string 
contribution to the CMB, for a given gauge group Ggut, depend on the 
superpotential coupling k, or equivalently on the symmetry breaking scale 
M associated with the inflaton mass scale, which coincides with the string 
mass scale. The relation between k and M is 



Mpi 2 TT ?/Q 



(80) 
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where 



dz 



z{{z + 1) ln(l + z-i) + {z-l) ln(l - z'^)}] . (81) 



The total quadrupole anisotropy has to be normahsed to the COBE data. 

A detailed study has been performed in Ref. [55] ■ It was shown that 
the cosmic string contribution is consistent with the CMB measurements, 
provided [82 

M ^ 2 X lO^^GeV ^ k<7x 10"^ . (82) 

In FiglH one can see the contribution of cosmic strings to the quadrupole 
anisotropy as a function of the superpotential coupling k |82j. The three 
curves correspond to A/" = 27 (curve with broken line), Af — 126 (full line) 
and J\f = 351 (curve with lines and dots). 




Fig. 4. Contribution of cosmic strings to the quadrupole anisotropy as a function of 
the superpotential coupling k. The three curves correspond to A/" = 27 (curve with 
broken line), TV = 126 (full line) and A/" — 351 (curve with lines and dots) [82) . 



The constraint on k given in Eq. (|82p is in agreement with the one found 
in Ref. [84]. Strictly speaking the above condition was found in the context 
of SO(IO) gauge group, but the conditions imposed in the context of other 
gauge groups are of the same order of magnitude since M is a slowly varying 
function of the dimensionality J\f of the representations to which the scalar 
components of the chiral Higgs superfields belong. 

The superpotential coupling k is also subject to the gravitino constraint 
which imposes an upper limit to the reheating temperature, to avoid gravitino 
overproduction. The reheating temperature Trh characterises the reheating 
process via which the Universe enters the high entropy radiation dominated 
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phase at the end of the inflationary era. Within the minimal supersymmetric 
standard model and assuming a see-saw mechanism to give rise to massive 
neutrinos, the reheating temperature is [82j 

Trh « ^^^{rM^xY'' , (83) 

where F is the decay width of the oscillating inflaton and Higgs fields into 
right-handed neutrinos 

1 /M,^' 



with TTiinfl — ^/2kM the inflaton mass and Ali the right-handed neutrino 
mass eigenvalue with Ali < m-infi/2. Equations (|5n|) . ((551) . (EH) lead to 

In order to have successful reheating, it is important not to create too many 
gravitinos, which imply the following constraint on the reheating tempera- 
ture [5S] Trh < lO^GeV. Since the two heaviest neutrinos are expected to 
have masses of the order of M3 ~ 10^^ GeV and M2 ~ 2.5 x lO^^ QgV re- 
spectively [HS], Mi is identified with Mi ~ 6 x lO'' GeV fSBj. The gravitino 
constraint on k reads [HI] k ^ 8 x 10"'^, which is clearly a weaker constraint 
than the one imposed from the CMB data. 

Concluding, F-term inflation leads generically to cosmic string formation 
at the end of the inflationary era. The cosmic strings formed are of the GUT 
scale. This class of models can be compatible with CMB measurements, pro- 
vided the superpotential coupling is smaller than 10^^ . This tuning of the 
free parameter k can be softened if one allows for the curvaton mechanism. 

According to the curvaton mechanism [88 89, , another scalar field, called 
the curvaton, could generate the initial density perturbations whereas the 
inflaton field is only responsible for the dynamics of the Universe. The cur- 
vaton is a scalar field, that is sub-dominant during the inflationary era as 
well as at the beginning of the radiation dominated era which follows the in- 
flationary phase. There is no correlation between the primordial fluctuations 
of the inflaton and curvaton fields. Clearly, within supersymmetric theories 
such scalar fields are expected to exist. In addition, embedded strings, if they 
accompany the formation of cosmic strings, they may offer a natural curva- 
ton candidate, provided the decay product of embedded strings gives rise to 



^ The linear mass density fj, gets a correction due to deviations from the Bogo- 
mol'nyi limit, which may enlare [87] the parameter space for F-term inflation. 
Note that this does not hold for D-term inflation, since then strings are BPS 
(Bogomornyi-Prasad-Sommerfield) states. 



30 Mairi Sakellariadou 



Contribution (%) ^^^q n=126 
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Fig. 5. The cosmic strings (dark gray), curvaton (liglit gray) and inflaton (gray) 
contributions to tlie CMB temperature anisotropies as a function of tlie tlie 
initial value of the curvaton field ■(/'init, and the superpotential coupling n, for 
N = 126 i82J. 

a scalar field before the onset of inflation. Considering the curvaton scenario, 
the coupling k is only constrained by the gravitino limit. More precisely, as- 
suming the existence of a curvaton field, there is an additional contribution to 
the temperature anisotropies. The WMAP CMB measurements impose [52] 
the following limit on the initial value of the curvaton field 

V-init ^ 5 X 10^3 (^) GeV , (86) 

provided the parameter n is in the range [10~^, 1] (see. Fig. [5]). 

The above results hold also if one includes supergravity corrections. This 
is expected since the value of the inflaton field is several orders of magnitude 
below the Planck scale. 
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D-term Inflation The early history of the Universe at energies below the 
Planck scale is described by an effective N=l supergravity (SUGRA) theory. 
Inflation should have taken place at an energy scale V^^"^ ^ 4 x 10^^ GeV, 
implying that inflationary models should be constructed in the framework of 
SUGRA. 

However, it is difficult to implement slow-roll inflation within SUGRA. 
The positive false vacuum of the inflaton field breaks spontaneously global 
supersymmetry, which gets restored after the end of inflation. In supergravity 
theories, the supersymmetry breaking is transmitted to all fields by gravity, 
and thus any scalar field, including the inflaton, gets an effective mass of the 
order of the the expansion rate H during inflation. 

This problem, known as the problem of Hubble-induced mass, originates 
from F-term interactions — note that it is absent in the model we have de- 
scribed in the previous subsection — and thus it is resolved if one considers 
the vacuum energy as being dominated by non-zero D-terms of some super- 
fields [90|91j . This result led to a dramatic interest in D-tcrm inflation, since 
in addition, it can be easily implemented within string theory. 

D-term inflation is derived from the superpotential 

Wi^fl = XS<P+-P- ; (87) 

S,'P^,'P+ are three chiral superfields and A is the superpotential coupling. 
D-term inflation requires the existence of a nonzero Fayet-Iliopoulos term 
^, which can be added to the Lagrangian only in the presence of an extra 
U(l) gauge symmetry, under which, the three chiral superfields have charges 
Qs = 0, Q^,^ = +1, and Q$_ = —1. This extra U(l) gauge symmetry can 
be of a different origin; hereafter we consider a nonanomalous U(l) gauge 
symmetry. Thus, D-term inflation requires a scheme, like 

GgutxU(1)^HxU(1)^^H^Gsm . (88) 

The symmetry breaking at the end of the inflationary phase implies that 
cosmic strings are always formed at the end of D-term hybrid inflation. To 
avoid cosmic strings, several mechanisms have been proposed which cither 
consider more complicated models or require additional ingredients. For ex- 
ample, one can add a nonrenormalisable term in the potential '92], or add 
an additional discrete symmetry [93j , or consider GUT models based on non- 
simple groups [51] , or introduce a new pair of charged superfields [HS] so that 
cosmic string formation is avoided at the end of D-term inflation. In what fol- 
lows, we show that standard D-term inflation followed unavoidably by cosmic 
string production is compatible with CMB data, because the cosmic string 
contribution to the CMB data is not constant nor dominant. Thus, one does 
not have to invoke some new physics. 
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In the global supersymmetric limit, Eqs. ([TS]). ([57)1 lead to the following 
expression for the scalar potential 



where g is the gauge coupling of the U(l) symmetry and ^ is a Fayet-Iliopoulos 
term, chosen to be positive. 

In D-term inflation, as opposed to F-term inflation, the inflaton mass ac- 
quires values of the order of Planck mass, and therefore, the correct analysis 
must be done in the framework of SUGRA. The SSB of SUSY in the infla- 
tionary valley introduces a splitting in the masses of the components of the 
chiral superfields 'P±. As a result, we obtain [96] two scalars with squared 
masses m\. = A^jS'p exp (jS'p/Mpj) ± and a Dirac fermion with squared 
mass TTif — A^jS'pexp (jS'p/Mpj). Calculating the radiative corrections, the 
effective scalar potential for minimal supergravity reads [82l96j 



As it was explicitely shown in Refs. '82 96 , D-term inflation can be com- 
patible with current CMB measurements; the cosmic strings contribution to 
the CMB is model-dependent. The results obtained in Refs. |82I96| can be 
summarised as follows: (i) g ^ 2 x 10~^ is incompatible with the allowed cos- 
mic string contribution to the WMAP measurements; (ii) for g ^ 2 x 10^^ the 
constraint on the superpotential coupling A reads A ^ 3 x 10^^; (iii) SUGRA 
corrections impose in addition a lower limit to A; (iv) the constraints induced 
on the couplings by the CMB measurements can be expressed as a single 
constraint on the Fayet-Iliopoulos term ^, namely ^ 2 x 10^^ GeV. They 
are shown in Fig. [B]. 

Assuming the existence of a curvaton field, the fine tuning on the couplings 
can be avoided provided [82l96j 



Clearly, for smaller values of A, the curvaton mechanism is not necessary. We 
show in Fig. [7] the three contributions as a function of ipi-ait, for A = 10~^ 
and g = 10^^. There are values of V'init which allow bigger values of the 
superpotential coupling A and of the gauge coupling than the upper bounds 
obtained in the absence of a curvaton field. 



(89) 
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Fig. 6. Cosmic string contribution to the CMB temperature anisotropies as a func- 
tion of the superpotential couphng A for different values of the gauge couphng g. 
The maximal contribution allowed by WM AP is represented by a dotted line [82l96j . 
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Fig. 7. The cosmic strings (dark gray), curvaton (light gray) and inflaton (gray) 
contributions to the CMB temperature anisotropies as a function of the the initial 
value of the curvaton field ijjinit, for A = 10~^ and g = lO"'^ [96] . 

Concluding, standard D-term inflation always leads to cosmic string for- 
mation at the end of the inflationary era; these cosmic strings are of the 
grand unification scale. This class of models is still compatible with CMB 
measurements, provided the couplings are small enough. As in the case of F- 
term inflation the fine tuning of the couplings can be softened provided one 
considers the curvaton mechanism. In this case, the imposed CMB constraint 
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on the initial value of the curvaton field reads [82196] 

V^init ^ 3 X f (^) GeV , for A E [10-\ IQ-^] . (92) 

The above conclusions are still valid in the revised version of D-term 
inflation, in the framework of SUGRA with constant Fayet-Iliopoulos terms. 
In the context of N=l, 3+1 SUGRA, the presence of constant Fayet-Iliopoulos 
terms shows up in covariant derivatives of all fermions. In addition, since the 
relevant local U(l) symmetry is a gauged R-symmetry |97], the constant 
Fayet-Iliopoulos terms also show up in the supersymmetry transformation 
laws. In Ref. ^8] there were presented all corrections of order g^/Mpj to the 
classical SUGRA action required by local supersymmetry. Under U(l) gauge 
transformations in the directions in which there are constant Fayet-Iliopoulos 
terms ^, the superpotential must transform as [97] 

= -^{^W , (93) 

otherwise the constant Fayet-Iliopoulos term ^ vanishes. This requirement is 
consistent with the fact that in the gauge theory at Mp\ oo the poten- 
tial is U(l) invariant. To promote the simple SUSY D-term inflation model, 
Eq. ([87]) . to SUGRA with constant Fayet-Iliopoulos terms, one has to change 
the charge assignments for the chiral superfields, so that the superpoten- 
tial transforms under local R-symmetry [55] . In SUSY, the D-term potential 
is neutral under U(l) symmetry, while in SUGRA the total charge of <P± 
fields does not vanish but is equal to — ^/Mpj. More precisely, the D-term 
contribution to the scalar potential V [see Eq. ([89]) ] . should be replaced by 
(gV2)(g+|<^+P+q_|0_|2+e)2 where 

<Z±-±l-p±-|^ with p++p^ = l. (94) 

In addition, the squared masses of the scalar components (/)± become 

ml = A2|5|2exp(|5|VMpi) ±3'e9± ; (95) 

the Dirac fermion mass remains unchanged. 

For the limits we imposed on the Fayet-Iliopoulos term ^, the correction 
^/Mpj is ~ 10~^, implying that the constraints we obtained on g and A, or 
equivalently on y^, as well as the constraint on ■0init still hold in the revised 
version of D-term inflation within SUGRA [63] . 

It is important to generalise the above study in the case of nonmini- 
mal SUGRA [53], in order to know whether qualitatively the above picture 
remains valid. A recent study [99] has shown that non- minimal Kahler poten- 
tial do not avoid the fine tuning, since the cosmic string contribution remains 
dominant unless the couplings and mass scales are small. For example, taking 
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into account higher order corrections in the Kahler potential, or considering 
supergravity with shift symmetry, we have obtained [99] that the 9% con- 
straint in the allowed contribution of cosmic strings in the CMB spectrum 
implies 

,/! < 2 X lO^^GeV <^ < 8 X 10"'^ . (96) 

The cosmic string problem can be definitely cured if one considers more com- 
plicated models, for example where strings become topologically unstable, 
namely semi-local strings. 

4 Cosmic Superstrings 

At first, and for many years, cosmic strings and superstrings were considered 
as two well separated issues. The main reason for this clear distinction may 
be considered the Planckian tension of superstrings. If the string mass scale 
is of the order of the Planck mass, then the four-dimensional F- and D-string 
self gravity is G4^f ~ ^idl) and G4/XD1 ~ C(ffs) (ffs stands for the string 
coupling), respectively, while current CMB measurements impose an upper 
limit on the self gravity of strings of G/i < 10~^. Moreover, heavy superstrings 
could have only been produced before inflation, and therefore diluted. In 
addition, Witten showed [100] that, in the context of the heterotic theory, 
long fundamental BPS strings are unstable, thus they would not survive on 
cosmic time scales; non-BPS strings were also believed to be unstable. 

At present, the picture has been dramatically changed (for a review, see 
e.g., Ref. |101j ). In the framework of braneworld cosmology, our Universe 
represents a three-dimensional Dirichlet brane (D3-brane) on which open 
fundamental strings (F-strings) end |102| . Such a D3-brane is embedded in 
a higher dimensional space, the bulk. Brane interactions can unwind and 
evaporate higher dimensional branes, leaving behind D3-branes embedded 
in a higher dimensional bulk; one of these D3-branes could play the role 
of our Universe [103]. Large extra dimensions can be employed to address 
the hierarchy problem 104j , a result which lead to an increasing interest 
in braneworld scenarios. As it has been argued [1051106] D-brane-antibrane 
inflation leads to the production of lower-dimensional D-branes, that are 
one-dimensional (D-strings) in the noncompact directions. The production of 
zero- and two-dimensional defects (monopoles and domain walls, respectively) 
is suppressed. The large compact dimensions and the large warp factors can 
allow for superstrings of much lower tensions, in the range between 10~^^ < 
G/i < 10^^. Depending on the model of string theory inflation, one can 
identify 107J D-strings, F-strings, bound states of p fundamental strings and 
q D-strings for relatively prime (p, q) , or no strings at all. 

The probability that two colliding superstrings reconnect can be much less 
than one. Thus, a reconnection probability P < 1 is one of the distinguish- 
ing features of superstrings. D-strings can miss each other in the compact 
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dimension, leading to a smaller V , while for F-strings the scattering has to 
be calculated quantum mechanically, since these are quantum mechanical 
objects. 

The collisions between all possible pairs of superstrings have been studied 
in string perturbation theory [108) . For F-strings, the reconnection proba- 
bility is of the order of g^. For F-F string collisions, it was found |il08j that 
the reconnection probability V is lO^'^ ^ 7^ 1. For D-D string collisions, 
10~^ ^ P 1. Finally, for F-D string collisions, the reconnection probability 
can take any value between and 1 . These results have been confirmed [109j 
by a quantum calculation of the reconnection probability for colliding D- 
strings. Similarly, the string self-intersection probability is reduced. When D- 
and F-strings meet they can form a three-string junction, with a composite 
DF-string. In IIB string theory, they may be found bound (p, q) states of p 
F-strings and q D-strings, where p and q are coprime. This leads to the ques- 
tion of whether there are frozen networks dominating the matter content of 
the Universe, or whether scaling solutions can be achieved. 

The evolution of cosmic superstring networks has been addressed numer- 
ically |110lllllll2lll3lll4lll5j and analytically pT6] . 

The first numerical approach [llOj . studies independent stochastic net- 
works of D- and F-strings, evolving in a flat spacetime. One can either evolve 
strings in a higher dimensional space keeping the reconnection probability 
equal to 1, or evolve them in a three-dimensional space with <C 1. These 
two approaches lead to results which are equivalent qualitatively, as it has 
been shown in Ref. [110] . These numerical simulations have shown that the 
characteristic length scale ^, giving the typical distance between the nearest 
string segments and the typical curvature of strings, grows linearly with time 



the slope C depends on the reconnection probability V , and on the energy 
of the smallest allowed loops (i.e., the energy cutoff). For reconnection (or 
intercommuting) probabihty in the range 10""^ ^ V ^ 0.3, it was found [llOj 



in agreement with older results [iP] . 

One can find in the literature statements claiming that ^(t) should be 
proportional to Vt instead. If this were correct, then the energy density of 
cosmic superstrings of a given tension could be considerably higher than that 
of their field theory analogues (cosmic strings). In Ref. 117J it is claimed 
that the energy density of longs strings pi evolves as pi — 2(a/a)p\ -~'P{p\/S,), 
where H = a/a is the Hubble constant. Substituting the ansatz ^(t) = 7(t)t, 
the authors of Ref. ^11 7J obtain 7 = — [l/(2i)](7 — V), during the radiation- 
dominated era. This equation has a stable fixed point at j(t) ^ V , implying 
that [117j ^ ~ Vt. However, Ref. [117j misses out the fact that intersections 
between two long strings is not the most efficient mechanism for energy loss 



(97) 



( (X ^ ^{t) cx ^/Vt , 
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of the string network. The possible string intersections can be divided into 
three possible cases (see, Fig. [Ij: (i) two long strings collide in one point 
and exchange partners with intercommuting probability Vi, (ii) two strings 
collide in two points and exchange partners chopping oflF a small loop with 
intercommuting probability V^] and (iii) one long string self- intersects in 
one point and chops off a loop with intercommuting probability V2, which 
in general is different than Vi- Only cases (ii) and (iii) lead to a closed 
loop formation and therefore remove energy from the long string network. 
Between cases (ii) and (iii) , only case (iii) is an efficient way of forming loops 
and therefore dissipating energy: case (iii) is more frequent than case (ii), and 
case (ii) has in general a smaller probability, since Vi ~ ■p2 |110j . However, 
the heuristic argument employed in Ref. ^117] does not refer to self-string 
intersections (i.e, case (iii)); it only applies to intersections between two long 
strings, which depend on the string velocity. However self-string intersections 
should not depend on how fast the string moves, a string can intersect itself 
even if it does not move but it just oscillates locally. 

The findings of Ref. [llOj cleared the misconception about the behaviour 
of the scale ^, and shown that the cosmic superstring energy density may be 
higher than the field theory case, but at most only by one order of magnitude. 

An important question to be addressed is whether cosmic superstrings 
can survive for a long time and eventually dominate the energy density of 
the Universe. This could lead to an overdense Universe with catastrophic cos- 
mological consequences. If the reconnection probability is too low, or equiva- 
lently, the strings move in a higher dimensional space and therefore miss each 
other even if V is high, then one may fear that the string network does not 
reach a scaling regime. The string energy density redshifts as where a 

stands for the scale factor. Since for a string network with correlation length 
^, there is about 1 long string per horizon volume, the string energy density 
is ~ yuja} . String interactions leading to loop formation guarantee a scaling 
regime, in the sense that strings remain a constant fraction of the energy 
density of the Universe. Loops do not feel the expansion of the Universe, so 
they are not conformally stretched and they redshift as Xjc?. As the loops 
oscillate, they lose their energy and they eventually collapse. Clearly, scaling 
is not a trivial issue for cosmic superstrings. 

In the first numerical approach |110j . where they have been only consid- 
ered independent stochastic networks of either F- or D-strings, it was shown 
that each such network reaches a scaling regime. This has been shown by 
either evolving strings in a higher dimensional space with intercommuting 
probability equal to I, or evolving strings in a three-dimensional space with 
intercommuting probability much smaller than 1. 

In a realistic case however, (p, q) strings come in very large number of dif- 
ferent types, while a (p, q) string can decay to a loop only if it self-intersects 
of collide with another (p, q) or (— p, —q) string. A collision between (p, q) 
and (p', g') strings will lead to a new string (p ± p', g ± g'), provided the end 
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points of the initial two strings arc not attached to other three-string vertices, 
thus they are not a part of a web. If the coUision between two strings can 
lead to the formation of one new string, on a timescale much shorter than 
the typical collision timescale, then the creation of a web may be avoided, 
and the resulting network is composed by strings which are on the aver- 
age nonintersecting. Then one can imagine the following configuration: A 
string network, composed by different types of (p, q) strings undergoes col- 
lisions and self-intersections. Energy considerations imply the production of 
lighter daughter strings, leading eventually to one of the following strings: 
(±1, 0), (0, ±1), ±(1, 1), ±(1, — 1). These ones may then self-intersect, form 
loops and scale individually. Provided the relative contribution of each of 
these strings to the energy density of the Universe is small enough, the Uni- 
verse will not be overdosed. 

This result has been confirmed by studying numerically the behavior of 
a network of interacting Dirichlet-fundamental strings {p, q) in Ref . |112j . To 
model {p, q) strings arising from compactifications of type IIB string theory, 
the authors studied [112] the evolution of nonabelian string networks. The 
positive element of such nonabelian networks is that they contain multiple 
vertices where many different types of string join together. Such networks 
have the potential of leading to a string dominated Universe due to tangled 
networks of interacting {p, q) strings that freeze. It was shown jll2j that 
such freezing does not take place and the network reaches a scaling limit. 
In this field theory approach however strings are not allowed to have differ- 
ent tensions, which is a characteristic property of cosmic superstrings. This 
issue has been addressed later in the context of modelling {p, q) cosmic su- 
perstrings |113j . It was found that such networks rapidly approach a stable 
scaling solution, where once scaling is reached, only a small number of the 
lowest tension states is populated substantially. An interesting question is to 
find out whether the field theory approach of Ref. [112] mimics the results of 
the modelling approach of Ref. |113j . Finally, performing full classical field 
theory simulations for a model of a string network with junctions, where the 
junctions can be thought of as global monopoles connected by global strings, 
it was shown ill5j that the evolution is consistent with a late-time scaling 
regime. Thus, the presence of junctions is not itself inconsistent with scaling. 

The cosmic superstring network is characterised [110' by two components: 
there are a few long strings with a scale-invariant evolution; the characteristic 
curvature radius of long strings, as well as the typical separation between two 
long strings are both comparable to the horizon size, (^{t) ~ \/Vt, and there 
is a large number of small closed loops having sizes t. Assuming there 
are string interactions, the network of long strings will reach an asymptotic 
energy density, where the energy density in long strings is 
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Thus, the fraction of the total density in the form of strings in the radiation- 
dominated era reads 

^ = . (100) 

Ptotal 3 V 

Recent numerical investigations [114] of strings evolving in a matter- or 
radiation-dominated FLRW background claim a weaker power law for the 
dependence of the scaling string energy density. More precisely, in Ref. [114j 
it was found that for V <; 0.1, the function p{\/V) is approximately flat, 
while for V ^ 0.1, the function pil/V) is well-fitted by a power-law with 
exponent O.GIq'Jj- The behaviour of the string energy density as a function 
of V has an important impact for the observational consequences of cosmic 
superstring networks. 

Oscillating string loops loose energy by emitting graviton, dilaton and 
Ramond-Ramond (RR) fields. Accelerated cosmic strings are sources of gravi- 
tational radiation, in particular from the vicinity of the cusps where the string 
velocity approaches the speed of light. Similarly, cosmic superstrings emit 
gravity waves but since the intercommutation probability is less than unity, 
their network is denser with more cusps, resulting in an enhancement of the 
emitted gravitational radiation. As it was pointed out [118j . the gravitational 
wave bursts emitted from cusps of oscillating string or superstring loops could 
be detectable with the gravitational- wave interferometers LIGO/VIRGO and 
LISA. 

One can place constraints on the energy scale of cosmic strings from the 
observational bounds on dilaton decays [119] . Considering that the dilaton 
lifetime is in the range lO^s ^ r < tdcc, one can obtain an upper bound 
r] ^ ■p~i/3 <^ lO^^GeV |110j for the energy scale of cosmic superstrings, which 
determines the critical temperature for the transition leading to string forma- 
tion. A lower reconnection probability allows a higher energy scale of strings, 
at most by one order of magnitude. 



5 Conclusions 

A realistic cosmological scenario necessitates the input of high energy physics, 
implying that models describing the early stages of the evolution of the Uni- 
verse have their foundations in both general relativity as well as high energy 
physics. Comparing the predictions of such models against current astro- 
physical and cosmological data one concludes to either their acceptance or 
their rejection, while in the first case one can also fix the free parameters 
of the models. One of the most beautiful examples in this interplay between 
cosmology and high energy phyics is the case of cosmic strings. 

Cosmic strings are expected to be generically formed during the evolu- 
tion of the Universe, provided the general theoretical picture we have in 
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mind is correct. However, many independent studies concluded in the ro- 
bust statement that cosmic strings have a hmited role in the measured CMB 
temperature anisotropics. Knowing the upper bounds on the contribution of 
strings to the CMB, one has to examine whether the theoretical models can 
be adjusted so that there is an agreement between predictions and data. This 
issue has been addressed in length here. In this respect, cosmology uses high 
energy physics to build a natural and successful cosmological model, while it 
offers back some means for testing high energy physics itself. 

Cosmic strings are a robust prediction of GUTs, or even M-theory. Even 
though their role in explaining the origin of the observed large-scale structure 
is sub-dominant, their astrophysical and cosmological implications remain 
important. Cosmic strings are a small but by no means negligible contribution 
to any successful cosmological model. 
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